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In traditional thermodynamics the Carnot cycle yields the ideal performance bound of heat engines 
and refrigerators. We propose and analyze a minimal model of a heat machine that can play a similar 
role in quantum regimes. The minimal model consists of a single two-level system with periodically 
modulated energy splitting that is permanently, weakly, coupled to two spectrally-separated heat 
baths at different temperatures. The equation of motion allows to compute the stationary power and 
heat currents in the machine consistently with the second-law of thermodynamics. This dual-purpose 
machine can act as cither an engine or a refrigerator (heat pump) depending on the modulation 
rate. In both modes of operation the maximal Carnot efficiency is reached at zero power. We 
study the conditions for finite-time optimal performance for several variants of the model. Possible 
realizations of the model are discussed. 
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I. INTRODUCTION 

Emerging technologies aim at operation based on quantum mechanics [1, 2], but power-supply or cooling devices are 
still governed by traditional (19th century) thermodynamics [3]. It is therefore imperative to examine the conceptual 
compatibility of the two disciplines as regards the performance of such devices [4-8]. On the practical side, progress 
in computer nanotechnology is currently constrained by the need to understand and optimize power and cooling 
generation on space-[9] and time- [10, 11] scales where quantum effects are unavoidable. This prompts the strive to 
achieve a better grasp on the fundamental thermodynamic bounds of quantum devices. 

In traditional thermodynamics, the Carnot cycle [3] yields the ideal performance bound of heat engines and pumps 
(refrigerators). In this cycle the evolution consists of "strokes" in which the system ("working fluid") alternates 
between coupling to the "hot" and "cold" heat baths. Yet in microscopic or nanoscopic devices, certainly when they 
operate quantum mechanically, such cycles pose a serious problem: on-off switching of system-bath interactions may 
strongly affect energy and entropy exchange, which casts doubts on the validity of commonly discussed models that 
ignore such effects. 

In this article we put forward a more rigorous approach to working cycles in quantum devices: we describe the 
steady-state dynamics of periodically- driven open quantum systems that are permanently coupled to heat baths by 
Floquet (harmonic) expansion of their coarse-grained Liouvillian evolution. The accuracy of this approach and its 
consistency with thermodynamics are ensured for weak system-bath coupling [12]. Here we apply this theoretical 
machinery to the description and performance analysis of a minimal model of a quantum heat machine (QHM), with 
the following features: 1) It is self-contained, i.e. described by a quantum mechanical Hamiltonian with the smallest 
possible number of degrees of freedom. 2) It is universal, i.e., it can act "on demand" as either a quantum heat engine 
(QHE) [6, 7, 13] that produces work, or as a quantum refrigerator (QR) [14, 15] that refrigerates a bath with finite 
heat-capacity, depending on a control parameter. 3) It is broadly adaptable to the available bath/environment or 
temperature. 

Our minimal model, the simplest to date, consists of a single qubit permanently attached to both baths and 
controlled by a harmonic-oscillator "piston" . A related (mostly numerical) study of QR has employed a qubit that 
alternately couples to one or another of two baths, whose spectra have different cutoffs [6] . Our goal is the development 
of a comprehensive analytical theory for universal, dual-purpose (QHE or QR) operation and its performance bounds 
in our minimal QHM model, wherein the spectral separation of the two baths plays a key role. Broad applicability is 
here ensured by bath engineering: attaching a "doorway mode" to an arbitrary bath acts as a bandpass filter that can 
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impose the required spectral separation on the two baths. Such doorway (filter)modes are realizable by interfacing 
the system (qubit) with the baths through a tunable cavity [16] or an impurity/defect in a periodic structure/chain 
[17]. 

The simplicity of the model allows closed-form analytical solutions, in which the piston-qubit coupling strength is the 
"knob" that can transform a QR to QHE (or vice versa) and controls their efficiency, after appropriately engineering 
the baths in question. Remarkably, Carnot efficiency is analytically shown to be achievable at the value of the control 
parameter that transforms the QR into a QHE. The optimal power and efficiency for finite-time cycles are analytically 
shown to surpass the established Curzon-Ahlborn bound [18]. 

In Sec. II the model and the analysis framework are introduced. In Sec. Ill the periodically modulated steady state 
is evaluated by the Floquet expansion method. In Sec. IV the steady-state thermodynamic relations are derived. In 
Sec. V we investigate the operation modes and bounds of the QHM, based on the steady-state solution of Sec. IV. 
In Sec. VI we investigate the finite-time optimal performance bounds. Realizations are discussed in Sec. VII. The 
conclusions are presented in Sec. VIII. 

II. MODEL AND TREATMENT PRINCIPLES 



The Hamiltonian of the QHM in question can be written as 

Hqhm = H s (t) + H B + H S b] (1) 
H SB = o- x (B H + B c ). (2) 

Here the control two-level system (TLS) is weakly coupled simultaneously to two baths via Hsb, where a x is the 
spinor x-component, B H and B c are respectively the operators of a very large hot bath (H) and of a finite cold bath 
(C). The TLS frequency is periodically modulated about its resonance frequency wo by the Hamiltonian 

Hs(t) = \o zV (t) (3) 

This model Hamiltonian is realizable by adiabatically eliminating a highly detuned level of a three-level system and 
allowing for a periodic AC Stark shift by a tio=me-dependent control field (Sec. VIIc). The fully quantized version 
of this model, wherein the classical time-dependent control field is replaced by a quantum harmonic-oscillator field 
dispersively coupled to the TLS, merits separate discussion. 
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A scenario that illustrates the model (Fig. 1-inset) is as follows: A charged quantum oscillator in a double-well 
potential which is "sandwiched" between the baths, a C-bath with finite heat-capacity and a nearly-infinite H-bath 
which serves as heat dump. The oscillation is periodically modulated, e.g. by off-resonant 7r-pulses, . These phase-flips 
control the heat current between the baths via the particle. 

In the refrigerator mode, corresponding to C-bath cooling, this model is reminiscent of the so-called sideband 
cooling: an optical Raman process in solids and molecules [19, 20]. Here the red- and blue- shifted TLS frequencies 
play the role of Stokes and anti-Stokes lines of sideband cooling respectively: heat is pumped into an upshifted line 
in the H-bath spectrum, at the expense of a downshifted C-bath spectral line, the energy difference is supplied by the 
modulation. In the engine mode, the opposite occurs: the modulation converts part of the heat-flow energy from the 
H-bath to the C-bath into work extractable by the control field. This entails energy transfer from the H-bath to the 
field. 

Work-extraction or refrigeration conditions are determined by the direction of power and heat flow (the heat 
current). The heat current, in turn, is given by the polarization rate of the TLS, obtained from the steady-state 
solution of a master equation (ME) for the TLS density operator. This ME, which allows for non-Markovian (bath- 
memory) effects, is accurate to second order in the system-bath coupling, at any temperature, as verified by us both 
theoretically [21] and experimentally [22, 23]. Inaccuracies of the ME[24] are negligible for weak coupling. (Sec Suppl. 
Mat. [25]) 

The refrigeration of a finite-capacity C-bath represents a succession of tiny temperature changes over many modu- 
lation cycles. Hence, the Born approximation underlying the ME is consistent with such cooling. The finite-capacity 
C-bath is assumed to have a continuous spectrum, since bath-mode discreteness and the associated recurrences may 
render the Born approximation invalid and preclude bath thermalization [26]. 

The TLS density matrix ps is assumed diagonal in the energy-state basis: starting at equilibrium, where off- 
diagonal elements of ps are absent, they remain so, when the TLS is subject to az~ modulation in Eq. (2). The 
interlevel transition rates and their non-Markovian time-dependence embody the quantumness of the ME. Sufficently 
fast periodic modulation of the TLS frequency at intervals r can change the detailed balance of the transition rates 
(quanta absorption and emission) and thereby allows either to heat up or cool down the TLS depending on r [10, 23]. 
Here our goal is the analysis of work extraction or heat flow between the baths enabled by periodic modulation of 
the TLS. Two alternative methods yield the same equations for these processes: a) Floquet (harmonic) expansion 
of the non-Markovian ME [10, 21] under temporal averaging (coarse-graining over a modulation period); b) Floquet 
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FIG. 1. Inset: Illustration of the QHM by double-well qubits (with periodically-modulated tunneling barrier) embedded 
between cold and hot baths, a) Sinusoidal-modulation effects (Eqs. (17), (18)). Harmonic (Floquet) peaks of the response 
superimposed on rising cold and hot bath spectra G 3 (uj) = A j lo 3 (top); and on the spectra of the same baths recalculated in 
the presence of different filter modes, transforming these baths spectra into skewed Lorentzians (according to Eqs. (39), (40)). 
The filtered spectra obey condition Al (left) or A2 (right). The parameters (in arbitrary units) used to calculate the spectra 
are: A c = 1, A H = 1/10, for the graph on the Mt:jf = 22, 7^ = 1, uf = 13 and uf = 1 and on the right jf = 1, jf = 2, 
uf — 13 and uif = 20. b) Same, under condition B (phase-flip modulation), for rising hot-bath spectra and cold-bath spectrum 
with cutoff. 
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FIG. 2. a) Currents and power as function of the machine modulation for an ideal heat machine. For A < A cr it operates as a 
heat engine and for A > A cr as a refrigerator, b) Efficiency for the same machine. The Carnot bound is reached at A = A cr . 
Inset: Periodic modulation of 3-level impurities embedded between hot and cold baths. The modulating field is detuned from 
the |e) — \u) transition by A u and exerts periodic AC Stark shift on |e). 

expansion of the Markovian evolution superoperator [12, 27]. 

III. FLOQUET EXPANSION AT STEADY STATE 
A. Non-Markovian master-equation approach 

Since ps is diagonal (see above) in the energy basis of the TLS (|e), \g)), the diagonal Hs(t) (Eq. (3)) yields the 
following rate equations from the non-Markovian ME [21] 

Pee(t) = -Pgg(t) = Rg{t)p g g " R e (t)p ee , (4) 

The non-Markovian, time-dependent |e) — > \g) and \g) — > |e) transition rates are given , respectively, by the real 
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part of the integrals [21, 28] 

R e (t) = 2Rc [ dt' exp[iw (t - t')]e{t)e* (t')<f>(t - t'), 
Jo 

R g (t) = 2 Re / dt' exp[-iuj (t - t')]e* '(t)e(t')<f>(t - t'), (5) 
Jo 

where e(t) is the periodically modulated phase factor (a unimodular periodic complex function), uj is the TLS 
resonance frequency and the bath-response (autocorrelation) function <&(£) = J dwG(w) cxp(— iivt) is the Fourier 
transform of to the bath coupling spectrum Gt(w). 

The TLS evolution caused by the baths under weak coupling conditions is much slower than their memory (cor- 
relation) time t c . Hence, in steady state, we can use time-averaged (coarse-grained) level populations and transition 
rates (See Suppl. Mat. [25]. In Eq.(3) the coarse-grained dynamics yields the following additive contributions of the 
two baths (labeled by j — C, H) to the harmonic expansion (labeled by m) of the time-averaged TLS polarization 



Cj — Pac—Pgg 



rfj(m) _ jjj(m) 

-(Ri^+R^)S+ H9 2 He 
where the time-averaged transition rates are found from the Floquet expansion of the modulation v{t) to be[3, 21, 28] 



R 3 ^ = 2nJ2 PmG j [±(uj + mA)]; 

m 

(6) 

Here 

Prn = \e m \ 2 , S m = - [ T e itiW)-«o)dt> e imAt dt (?) 

T Jo 

are the probabilities of shifting the j-th-bath coupling spectrum G^co) by mA, A = from the average frequency 
loo, and 

/+oo 
e^iBj^BjiO^dt = e w / T G J (-w) (8) 
-OO 



For a bosonic bath (h = ks = 1):: 



G>) =GS(w)(n j (w) + 1); G» = | 5 »| 2 p»; (9) 

gj(oj) being the system-bath coupling, pj{w) the bath -mode density and nj(cv) = — ) the w-mode thermal 
occupancy. 
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B. Markovian master equation approach 

An alternative method is based on the combination of the weak-coupling limit and Floquet expansion [27] of 
a periodically-flipped qubit coupled to two baths is based on the Lindblad-Gorini-Kossakowski-Sudarshan (LGKS) 
operator, as expounded in the tutorial in Rcf. [12]. It has the advantage of ensuring positivity and additivity of the 
evolution due to the two baths. The Markovian master equation in the interaction picture reads: = Cp where Cj. m 
we expand C = ^,£ J m , where again j — H, C, and m is the Floquet harmonic. The expansion yields 



o3(ff( W0 + mA)([O,a + ] + [<7-,p<7+]) + 

G'(-Lu -mA)([a + p,<j-} + [<7+ pa"])) (10) 

This (LGKS) approach yields the same equations of motion for the qubit polarization as the time-averaged non- 
Markovian approach above (Eq. (5)-(6)). 



IV. STEADY-STATE THERMODYNAMIC RELATIONS 



The qubit steady-state (found from Eq.(10)) has the diagonal form characterized by the population ratio w 



P = 



/ P~ee 



V Pgg J 



w = 



^ qJ P m &(uJo + mA)e~ 



Pgg 



(11) 



where, as in Eq.(6), A = — . The cold (hot) current is then given by 



J c(H) = E( w o + ™A)P m G c W(c + mA) 6 T ° {H '~ W (12 ) 
* — ' w + 1 

m 

The magnitudes and signs of these steady-state currents are the same as for the time-averaged (non-Markovian) 
ME solutions detailed above (Eq. (5)-(6)). 

The first law of thermodynamics allows to define the stationary power as 

V = -(J c + Jh). (13) 



i.e. the power investment V by the piston is negative when the machine acts as an engine. One obtains from (13) 

(14) 



p = £("o + mA)P„ 



w + l 



G : >(uj + mA)(w-e T o ) 



One can show using the Spohn theorem for Markovian evolution [29] that the heat currents satisfy the second law 
of thermodynamics 



S{t) = -Tr(p(t) Inp(t)) , j t S(t) ^ + ^ > (15) 



where the left hand side of the inequality (15) is the entropy production rate. 

These standard thermodynamical relations imply the validity of the Carnot bound on the engine efficiency rj, and 
on the coefficient of performance (COP) for the refrigerator 

-V . 1 T c 

V= -j- < I" 7p~, 

COP=^-<-^—. (16) 

V. UNIVERSAL MACHINE OPERATION MODES 

The choice of parameters that may affect the QHM operation are mainly the shape of the modulation v(t) and 
the form of bath-response (coupling spectra) spectral densities G c {lo) and G H (w). These choices may enable the 
machine to act as both an engine and a refrigerator, as shown below. In the following we discuss two such choices of 
the modulation and the requirements each type entails on the bath spectra. 

A. Sinusoidal modulation 

We consider the sinusoidal time-dependence of the external (modulating) held, i.e. 

uj(t) = w + AA sin(Ai) (17) 

under the condition 

0<A«1. (18) 
The condition (18) implies that only the harmonics m = 0, ±1 with 

P m =o - 1 - y - p m=±i - y (19) 

should be taken into account. Even under this simplifying condition the formulae for heat currents and power arc 
complicated (see Appendix A). More detailed analysis of the various terms in those formulae shows that in order to 
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reach the Carnot bound, we have to reduce the number of relevant harmonics to two. This can be done by system - 
bath coupling engineerings (Sec. VII), so as to impose the spectral separation of the two baths, in the cases discussed 
below. 

Al We assume that the upper cutoff of C nearly coincides with lower cutoff of H and lo is near the two cutoffs (Fig. 
la): 

G c (lo) ~ for to > lu q , G h {lo) ~ for to < to a , (20) 

We discuss this case in more detail, as the condition (20) is the easiest to implement in practice (see Section VII). This 
condition leads to the following simplified formulae for heat currents and power obtained from the general expressions 
(Appendix A) 

J H = {to + A)^( e - { — } - e- { — } ), 
J c = -(wo - A)^(e~ { ^T' - e- { ^r>), 

V = -2AJf{e ( t h > - e ( T c >), (21) 

where the positive normalizing constant is 

, = ^ G c (co ~A)G H (io + A) f22) 

4 G c (co - A) [1 + e- { ^T ] ] + G H (co + A) [1 + e "(^)~ ' 
It follows from (21) that there exists a critical value of the modulation frequency 



^t^TFc' (23) 



such that for A < A cr the machine acts as an engine with the efficiency 



and for A > A cr as a refrigerator with 



2A 



C0 P=^-^. (25) 
2A v ; 



At A = A cr the engine/refrigerator reaches its maximal Carnot efficiency/COP. This corresponds to the vanishing 
value of power/cold current (see Fig. 2a). The operation mode change prevents the machine from breaking the second 
law. Similar behavior was obtained [30] by numerical calculation for a different model of a quantum machine. 
A 2 We assume the following conditions of the bath spectra (Fig. 1 a) 

G c (lj) ~ for to w lo ± A , G H (to) ~ for lo < lo . (26) 
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This condition yields (App. A) 



J H = {lo + A)^'{e~ ( ^r ) - e" ( ^ } ) 
J c = uo^'ie-^ - e~ ( ^r ] ) 

V = -Aj^'{e~ { ^r ) -e~ { ^>), (27) 



where the positive normalizing constant for this case is 



jyi _ ^_ G c (u Ja )G H (u Ja + A) g 

4 GC( Wo )[l + e - ( ^ ) ] +^G H (Lu + A)[l + e-^ ) ] 



From Eq. (27) it follows that the critical modulation frequency is 



A cr =Lu TH Tc . (29) 



Namely for A < A cr the machine acts as an engine with the efficiency 



and for A > A cr as a refrigerator with 



'-iS+A- (30) 



COP= "° a A . (31) 



At A — A cr the engine/refrigerator reaches its maximal Carnot efficiency/ COP. 
B 7r-flips modulation 

Periodic 7r-phase shifts (phase flips) with alternating sign give rise, to only two leading harmonics, corresponding 
to two symmetrically-opposite frequency shifts in the Floquet expansion of the probability distribution (since P = 
for symmetry reasons) [21, 28]. 

P±i « (2/tt) 2 ; 

G j (cj ) -^G j (co ±A) (32) 

If A, is comparable to l/t c , the inverse memory time of the cold bath, we may require that at lj ~ lu + A the TLS 
be coupled only to the H bath, while at co ~ cj — A it is coupled to both the C and the H bath. This is tantamount 
to the requirement that 
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G h (lj + A) » G c (co + A),G h (lj - A),G c (u - A) (33) 
This requirement can be satisfied if C has an upper cutoff 

Wcut < w + A. (34) 

For H, by contrast, G H is only required to rapidly rise with w, which is true for blackbody radiation in open space, 
G H (ui) oc uj 3 , or for phonons in bulk media. 

If Eq. (33) is satisfied we find the same steady-state expression for the currents and power and the same physical 
behavior as for case Al: A cr is given by Eq. (23), r\ by Eq. (24) and COP by Eq. (25). 

VI. FINITE-TIME OPTIMAL PERFORMANCE 

The vicinity of the critical frequency is not a useful working regime, since power and currents are small there. Much 
more important is the region of parameters where the power or cold current are maximal. We wish to find these 
parameters and the corresponding efficiency of our QHM, and compare it to the so-called Curzon-Ahlborn efficiency 
at maximum power of a macroscopic Carnot-type enginc[18] 

%m = 1-(|^) 1/2 (35) 

To this end we consider the case A2 (Eq. (26)) for our universal machine and compute the maximal power V which 
is produced for the optimal modulation frequency A max under the following simplifying assumptions. 

a) The spectral density G h (uj) is flat, i.e. -£jG H ~ near uj = uj + A max . 

b) high-temperature regime e 7 ^ + p 

Under these assumptions, the modulation A max that yields the maximal power and the corresponding efficiency 
Vmax are found to be 

_1 _(!-^) 
A max — -A cr , Tj max — — - ^ > tic a (36) 

i H 

In the Table 1 we present the results of cases Al (Eq. (20)) and A2 (Eq. (26)). While the relation A max — \ A cr is 
true for all cases, at maximal-power the efficiencies are different. In the case A2 the Curzon-Ahlborn bound is always 
exceeded. In fact, t\ca is then the minimum of (t]) max . Numerical calculation shows that the Curzon-Ahlborn bound 
is exceeded even if the temperatures are not high (A w wq/2, Th — 2,Tc — 1). 
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TABLE I. 



Case 




Efficiency (A < A cr ) 


Efficiency at maximum power 


Relation to Curzon-Ahlborn 


Al 


, .. T H -T C 
U °T H+ T C 


2A 


„ _ 2(T H -T C ) 
r/max — 3Th+ t c 


r)max < VCA 


A2 


^ Th t7 c 


A 

UQ + A 


„ _ T H -T C 
Umax — T H +T C 


Vmax > VCA 



VII. REALIZATION CONSIDERATIONS 

In the following we lay out general guidelines for experimental realizations i.e. the key considerations for satisfying 
the operating conditions discussed above: 

a) Spectral separation of coupling to the C and H baths via Debye cutoffs 

Case B (periodic 7r-flips) is the most flexible so far as spectral separation is concerned: it allows for an arbitrary 
rising coupling spectrum of the H bath, including a bulk-solid phonon bath or the blackbody radiation spectrum 
G h (uj) oc lo 3 . The C bath should have G c (cj = 0) 7^ 0, as for 1/f noise spectra [31]. Preferably, the C bath is to have 
an upper cutoff, as is the case for phonon baths in crystals (Debye cutoff) and Ohmic or super-Ohmic noise baths [31]. 
By contrast, a single mode of a cavity with finite linewidth (finesse) is a Lorentzian bath[32], whose lack of cutoff 
would lower the machine efficiency or COP, unless the upshifted (anti-Stokes) frequency uj + A is at the far tail of 
the Lorentzian (Fig. 1). 

The spectral separation requirements are stricter in cases Al and A2 (under sinusoidal modulation). To impose 
these requirements, we can choose the materials C and H in such a way that the Debye frequencies match the desired 
temperatures of the baths 

lj% ~ T c , uj% ~ T H . (37) 
The coupling spectra of the baths are assumed to have the Debye shape 

^M = /i(^) 8 rr^75j^-H) ( 38 ) 

where fj are bath-specific constants and 9 is the Heaviside step function. 

b) Filters 

In general, the satisfaction of the spectral separation conditions can be facilitated by imposing a "filter" onto the 
qubit-bath coupling spectrum. To this end, we consider the model whereby the qubit is (weakly) coupled to two 
harmonic-oscillator "filter" modes with resonance frequencies w/ each such mode in turn, is coupled to the respective 
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(C or H) bath via coupling spectrum G-?(w), i.e. each filter mode mediates between the qubit and the respective bath. 
The qubit is then effectively coupled to the filter-mode bath response via[17] 



/V 7T ( CJ _(^ /+ AiH)) 2 + (7TG^H) 2 ' 



where 7^ is the coupling rate of the qubit to the filter mode, and 



AiM=P(f°^'^9) (40) 

JO W — LO 

P being the principal value, is the respective bath-induced Lamb shift [17, 33]. The filter-mode response spectrum (39) 
is a "skewed Lorentzian" for a completely general spectrally structured G-?(w). In particular, a cutoff or a bandcdgc 
of G J (w) curtails the Lorentzian and makes it strongly skewed [17], whereas if G-?(w) is spectrally flat, it is a simple 
Lorentzian centered at ojj. Such "filtering" can suppress undesirable tails of G 3 (u) and thus enforce conditions (20) 
or (33) in a broader range of media and parameters (Fig. la). 

In addition, the desired bath spectra are achievable by engineering. We note recent advances in microcavities, 
photonic crystals and waveguides that can help reshape photon bath spectra, as well as their phonon-bath counterparts, 
e.g., periodic structures with acoustic bandgaps [13]. 

c) Qubit realization and modulation 

A qubit may be realized by atomic or molecular multilevel impurities embedded at the interface between different 
material layers (Fig. 2- inset) A qubit may also be realized by a symmetric double- well (DW) potential with two bound 
states |e), \g): for example, a single-electron quantum dot or a superconducting Josephson qubit [1, 2] as well as its 
ultracold-atom analog [34]. The symmetric and antisymmetric superpositions of |e) and \g), the eigenstates of a x (Eq. 
(2)), are localized on the left- and right-hand well, respectively. In either bath the Debye cutoff should conform to the 
temperature as specified above. One example is that of a DW quantum dot "sandwiched" between dielectric layers 
with different Debye cutoffs lo° d (Fig.l-inset) Another example is an ultracold-atom DW qubit embedded between two 
optical lattices, where phonons have different uP D (Fig. 2- inset). 

The required modulation that conforms to Eqs. (19) or (23) may be realized by changing the energy difference of 
|e) and \g) by time-dependent AC Stark shift[21]. If \u) is an upper state with energy uj u , then an off-resonant control 
field detuned by A u from the |e)-|u) resonance will realize the piston-system coupling in Eq. (3) with u(t) = Q 2 (t)/A u , 
Q(t) being the Rabi frequency of the control field. (Fig. 2-inset).Thus, in DW qubits realized by ultracold atoms, 



15 

the optically-induced potential barrier between the wells may be periodically modulated or flipped in sign [34] (Fig. 
1-inset) 

In order to illustrate the performance of such a machine, we numerically compute the dependence of the basic 
thermodynamical parameters on the rate A. As shown in Fig. 2, the operation mode, as well as its efficiency, depends 
on A. In the plotted example, the overlap of spectral densities of the two baths vanishes, i.e. they fully satisfy the 
spectral separation condition A, hence the ideal Carnot efficiency is reached. If the overlap does not vanish, the 
Carnot bound is not attained. 

VIII. CONCLUSIONS 

A TLS (qubit) with energy modulation has been shown to constitute a minimal model for a universal quantum 
heat (QHM). As shown, this machine can be switched at will from an engine mode to refrigerator mode and vice 
versa by varying the modulation rate. This could be useful in situations where both kinds of thermal machines are 
needed, and the operation is simplified by having a single machine. The rate also allows to control the efficiency of 
our machines and in this way keep them optimized under bath temperature change, which underscores the versatility 
of our machine. Recently we have shown that this machine may violate the unattainability of the absolute zero (the 
third law) in the refrigeration mode for certain bath models [25]. 

The switching from engine to refrigerator mode in this machine occurs at a critical rate, i.e. at the critical point 
where Carnot limit appears to be close to be broken, yet instead the mode switching prohibits breaking this limit. 
The maximum efficiency is reached at the critical point, where the machine yields zero power, consistently with the 
second law. Practical engines are however designed to yield maximum power. It is usually presumed that the upper 
bound is the Curzon-Ahlborn efficiency [18]. However, as our model shows, this bound can be exceeded, it is in fact 
the lower bound for our finite-time engine. 

Both impurities and double- well qubits embedded in appropriate environments may act, depending on their energy 
modulation, in cither the engine (QHE) or the refrigerator (QR) mode of the proposed QHM, with potentially 
significant technological advantages. In particular, an impurity or quantum-dot fast-modulated qubit "sandwiched" 
between two nanosize solid layers (Fig. 1-inset) may act as a nanoscopic refrigerator (heat-pump) of a transistor (chip), 
that is much more miniaturized and less power-consuming than currently available microelectronic refrigerators [9]. 
Under slower modulation rate, the same setup may act as electron-current generator without external voltage bias, 
and a substitute for phase-coherent electron-current control [35] (whereas in the present scheme the modulating field 
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need not be coherent). 

To conclude, the present scheme demonstrates the ability of systematic quantum analysis of driven open systems 
to yield accurate, physically lucid, expressions for steady-state heat machine performance. The analysis shows that 
quantum mechanics and thermodynamics can be fully compatible on the level of an elementary (single- qubit) system 
provided one correctly accounts (by the Floquet expansion) for its entropy and energy exchange with the baths. This 
essential point has not been properly accounted for by previous treatments of quantum heat machines. 
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Appendix A. Heat currents and power under sinusoidal modulation 



In general, under the conditions(17), (18), only three frequencies w = uj 0} ujo ± A need to be taken into account. 
The expressions for power , cold current and hot current (up to terms of the order A 2 read as follows. 

1) The power 



E^c^G'^o + mAJll + e *< ) 
(PiPo (g c (lu + A)G C M (e - ^ - e "*£) + G c (oj a + A)G h (oj ) (e - ^ - e _ *&) + 
G H (lj + A)G c (lu ) {eT^T - e"^) + G c (w - A)G c (uj q ) (e~^ - e~^) 
+G H (u - A)G C M (e~^ - e - ^) + G c {lu - A)G H '(wo) (e - ^ - e - ^) + 
G H (ui + A)G h (lj ) (e~^ - e"^) + G H \u> - A)G H (u Jo ) (eT% - e - ^)) + 
(g c (w + A)G C K - A) (e - ^ - e - ^) + G c \u + A)G h {oj - A) (e - ^ - e - ^) + 
G (cj + A)G (ijJq - A) (e~^T _ e -T^j + G (cj + A)G (uiq — A) {e~^T _ e -TjrJJJ 
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2) The hot current 

(pfaoG H (w )G c (lo ) (e"^ - e~^) + PxPq (-w G c (wo + A)G ff (wo) (e - ^ - e - ^) + 
(w + A)G H (w + A)G c (w ) (e - ^ - e - ^) 
-(w - A)G H (w - A)G c (w ) (e - ^ -e"^) + w G c (w - A)G ff (w ) (e~^ -e - ^) + 
AG H (w + A)G H (w ) (e~^T - + AG H (w - A)G H (w ) (e~*& - e"^)) + 

(-(w - A)G c (w + A)G H (w - A) (e~^ - e~^) + (w + A)G ff (w + A)G c (w - A) (e~^ - e"^) 

+2AG H (w + A)G H (w - A) (e~"**r - e~^) + 
(w - A)G H (w - A)G c (w - A) (e'^T - e - ^))) 

3) The cold current expression can be obtained from Jh by interchanging C and H. 
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